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This  is  one  step  in  a  program  aimed  at  classifying  solutions  oi  the 

'  \ 

Riemann  problem  for  2  x  2  .  hyperbolic  quadratic  conservation  laws.  \  Such 

A  A 

conservation  laws  approximate  a  general  2x2  system  of  conservation  laws  in 

/A  . s'  *-)C 

a  neighborhood  of  a  point  at  which  strict  hyperbolicity  fails.  We  give '.the 

r 

solution  for  the  symmetric  systems  in  Region  III  of  the  four  region 
classification  of  Schaeffer  and  Shearer.  The  solution  is  based  on  the 
qualitative  shape  of  the  integral  curves  described  by  Schaeffer  and  Shearer 
and  a  numerical  calculation  of  the  Hugoniot  loci  with  their  shock  types. 
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SIGNIFICANCE  AND  EXPLANATION 

A  2x2  system  of  conservation  laws  is  a  system  of  partial  differential 
equations  of  the  form 


( 1) 


lit  +  £(a)x  =  0 


where  u  =  (u,v),  f  =  f(g)  are  in  R  x  r  and  where  x  e  R,  t  >  0.  Such 

equations  arise  in  gas  dynamics,  elasticity,  oil  reservoir  simulation  and 

other  areas  of  engineering  when  diffusion  is  neglected.  In  solutions  of  ( 1 ) , 

information  travels  at  speeds  X.  and  X_  given  by  the  eigenvalues  of  the 
3f  2 

matrix  Since  this  matrix  depends  on  u,  v,  the  speeds 

rsj 

depend  on  the  solution,  and  this  leads  to  the  formation  of  discontinuities 

called  shock  waves.  The  present  paper  deals  with  the  classification  of 

Riemann  problem  solutions  (solutions  which  evolve  from  a  single  discontinuity 

3f 

at  time  t  =  0)  near  an  isolated  point  at  which  —  is  a 


X1  and  X2 
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THE  CLASSIFICATION  OF  SOLUTIONS  OF  QUADRATIC  RIEMANN  PROBLEMS  (II) 
E.  Isaacson*'6,7'8,9, 10  and  B.  Temple1 ' 2’  3' 4f 5 


§  1.  Introduction 

This  is  the  second  in  a  series  of  papers  [5]  in  which  we  give  the  solution  of  the 
Riemann  problem  for  quadratic  conservation  laws 


u  |  {a^2  +  2b ^v  +  c.jV2]^  »  0  , 

v  +  -  f a„u2  +  2b_uv  +  c.v2l  “  0  , 
t  2  2  2  2  x 

with  initial  data 


u(x,  0 ) 


(  Sl  5  (VV  ' 

\  a*  =-  (W  ' 


x  <  0 
x  >  0 


where  u  =  (u,v).  This  is  one  of  the  steps  in  a  program  outlined  in  [6], 

Solutions  of  such  conservation  laws  aproximate  the  solutions  of  a  general  2*2 
system  of  conservation  laws  in  a  neighborhood  of  an  isolated  point  at  which  strict 
hyperbolicity  fails. 

We  use  the  normal  form 

12  2 
ufc  +  j  {au  +  2buv  +  v  )x  *  0  » 

(  1) 

v^  +  j  {bu2  +  2uv(x  ■  0 
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of  Schaeffer  and  Shearer  [12] 


He  let  X^tu)  <  X2(u)  denote  the  eigenvalues  of  system  (1);  we  note  that  for  system 
(1),  X1  *  X2  only  when  u  ■  0,  in  which  case  X  ■  0. 

In  the  first  paper  [5],  we  presented  the  solution  of  the  Riemann  problem  for  the 
parameter  values 

a  >  2  ,  b  -  0  . 

This  corresponds  to  the  symmetric  systems  in  Region  IV  of  [12,5].  In  the  present  work  we 
give  the  solution  of  the  Riemann  problem  for  the  parameter  range 

1<a<2,  b  *■  0  , 

which  corresponds  to  the  symmetric  systems  in  Region  III.  We  refer  to  [S]  for  a  detailed 
discussion  of  the  problem  and  the  notation. 

We  obtain  the  shapes  and  shock  types  of  Hugoniot  loci  H(u_)  by  means  of  numerical 
calculation.  This,  together  with  the  qualitative  features  of  the  integral  curves  given  by 
[12,2]  (see  Fig.  1)  is  the  basis  for  our  construction  of  the  solutions.  We  present  the 
solution  of  the  Riemann  problem  in  a  series  of  diagrams  for  representative  values  of  u^ 
The  qualitative  features  of  the  solution  diagrams  change  precisely  when  u^  crosses 
H( 0)  or  crosses  a  ray  X1  ”  0.  In  Regions  I-III  the  Hugoniot  locus  H(0)  consists  of 
three  straight  lines  (called  axes),  while  in  Region  IV  this  locus  is  a  single  line  [12]. 
For  the  symmetric  systems  in  Regions  I-III  (i.e.,  a  <  2,  b  “  0),  one  of  the  axes  is  the 
u-axis  and  the  other  two  are  given  by 

v  “  ±  /2  -  a  u  . 

The  rays  X 1  *  0  are  given  by 

f  -  /a  u  in  the  lower  half  plane 
v  "  { 

+  /a  u  in  the  upper  half  plane  . 

we  give  the  solution  of  the  Riemann  problem  for  representative  values  of  uL  in  the  lower 

half  plane  (v  <  0)  because  of  the  symmetries  in  system  (1)  (see  [5]).  Since  the 

solutions  change  qualitatively  only  at  the  rays  H(0)  and  X^  -  0,  we  choose  a 

representative  value  of  u_  from  each  interior  and  bounding  ray  of  the  sectors  (see  Fig. 

~L 

2) 


-2- 


A1  ■ 

k 

ifl,2  <  9  <  0}  , 

A2  “ 

fu 

:  9.  <  0  <  912) 

A  3  “ 

fu 

:034  <  9  <  9.} 

A4  “ 

fa 

«-*  <  6  <  834} 

where 

912  «  arctan  ( — / 2  -  a)  , 

9  »  arctan  (+/2  -  a)  -  u  , 

34 

0#  «  arctan  (-/a)  . 

The  angles  correspond  to  the  axes  and  the  ray  “  0,  respectively. 

A  new  feature  in  Region  III  is  that  the  Hugoniot  locus  Htu^)  and  the  1-wave  curve 
W1 (uL)  can  be  disconnected  [5].  This  makes  the  solutions  more  complicated  than  those 
found  in  Region  IV. 

In  Section  2  we  describe  the  elementary  waves  and  Hugoniot  loci,  and  in  Section  3  we 
present  the  solution  of  the  Riemann  problem. 

We  wish  to  acknowledge  the  assistance  of  Dan  Marchesin  and  Bradley  Plohr  with  the 
numerical  aspects  of  the  problem.  We  also  thank  Michael  Shearer  and  David  Schaeffer  for 
helpful  discussions  on  many  aspects  of  quadratic  Riemann  problems.  Finally,  we  would  like 
to  thank  Dan  and  Paulo  paes-Leme  for  their  gracious  hospitality  during  our  stay  in  Brazil. 


§2.  Elementary  Waves 

The  rarefaction  waves  of  system  (1)  are  obtained  from  the  integral  curves  of  the 
eigenvector  fields  of  the  Jacobian  A(u) .  The  shock  waves  of  system  (1)  are  obtained  from 
the  Hugoniot  loci.  The  general  solution  of  the  Riemann  problem  that  we  construct  is 
obtained  by  composing  these  waves.  The  solution  consists  of  a  1 -composite  wave  followed 
by  a  2-composite  wave  where  a  p-composite  wave  is  a  succession  of  shock  and  rarefaction 
waves  of  the  p  family  [51.  The  p-wave  curve  U/ ^ ( u ^ >  consists  of  all  final  states  in  p- 
composite  waves  with  initial  state  u^. 

The  qualitative  shapes  of  the  integral  curves  for  the  symmetric  systems  (1)  in  Region 
III  are  depicted  in  Figure  1  (see  [12,2]). 

2. 1  Hugoniot  Loci 

The  structure  of  the  Hugoniot  loci  is  depicted  in  Figs.  3A-J  for  representative 
values  of  u^.  These  figures  indicate  the  general  topological  structure  of  the  loci  as 
well  as  the  location  of  shock  types.  As  uT  varies,  these  features  change  qualitatively 

~Li 

only  at  the  boundaries  of  the  Sectors  Aj  -  A4  (and  at  X2  ”  0).  The  topological 
structure  of  the  Hugoniot  locus  changes  at  the  axes,  while  at  the  lines  X  •  0  the 
topological  structure  remains  fixed  but  the  location  of  shock  types  changes.  (He  do  not 
include  X2  »  0  as  a  boundary  since  the  change  in  shock  types  at  Xj  •  0  does  not  affect 
the  Riemann  problem  solutions.  Specifically,  the  wave  curves,  and  hence  solutions,  do  not 
change  qualitatively  at  Xj  »  0. ) 

Note  that  the  Hugoniot  loci  of  states  in  quadrant  IV  of  the  u,v-plane  determine  all 
Hugoniot  loci  because  the  reflection  of  a  Hugoniot  locus  about  either  the  u-  or  v-axis  is 
again  a  Hugoniot  locus  (5).  For  clarity,  we  include  the  Hugoniot  loci  for  sector  and 

its  boundary  rays  9  «  9  -  9J4  because  these  are  relevant  to  the  solution  of  the 


Riemann  problem. 


3.  Solution  of  the  Rlemann  Problem 

In  Figs.  4A-J  we  give  the  solution  of  the  Riemann  problem  for  the  symmetric  systems 
(1)  In  Region  III.  In  each  diagram,  Is  fixed  and  an  arbitrary  point  represents  Uj^. 

Figs.  4C,  E,  G  and  I  depict  the  solutions  for  uL  in  Sectors  A^,  Aj»  A3  and  A4* 
respectively,  in  each  of  which  the  solution  diagrams  are  qualitatively  the  same.  For 
completeness,  we  include  solution  diagrams  for  u.  on  the  boundaries  of  these  sectors: 

The  ray  9  *  S12  is  the  axis  separating  Sectors  A3  and  A2»  and  the  ray  9  “  93<  is 
the  axis  separating  Sectors  A3  and  A^  The  ray  9  «  9.  which  separates  Sectors  A2 
and  A3  is  the  key  ra;,  on  which  “  0. 

In  each  solution  diagram,  the  1-wave  curve  (^(u^)  consists  of  the  union  of  the  1- 

shock,  1-rarefaction,  and  1-composite  curves  for  u^.  (See  the  Legend.)  For  each  point 

u  e  W.(uT  ),  tve  portion  of  the  2-wave  curve  W.,(u_)  appearing  in  the  diagrams  is 
~m  1  £  mb 

obtained  by  proceeding  from  u  along  2-ahock  and  2-rarefaction  curves  as  far  as  possible 
in  the  direction  of  the  arrows. 

The  solution  of  the  Riemann  problem  <  u, ,u„  >  consists  of  a  1-wave  with  left  state 

~K 

u  and  right  state  u  followed  by  a  2-wave  with  left  state  u_  and  right  state 

where  the  intermediate  state  u  is  determined  as  follows:  start  from  u_  and  follow 

~rn  *wk 

the  2-wave  curve  backwards  from  uR  (opposite  the  direction  of  the  arrows)  until  you 
reach  a  point  u  in  W,(u,  )•  The  state  u  so  constructed  satisfies  u_  e  W_(u_)  and 
defines  the  waves  in  the  solution.  This  procedure  is  not  well-defined  precisely  when  uR 
is  in  either  the  compressive  portion  of  H(u^)  or  the  triple  shock  curve.  (The  triple 
shock  curve  is  [ABjJ  in  Figs.  4C-E  and  is  [-u#A]  in  Figs.  4F,  G.  There  is  no  triple 
shock  curve  in  the  other  figures.)  In  the  case  of  ambiguity,  the  solutions  are  unique  ir 
the  x,t-plane  since  then  all  shock  speeds  that  occur  are  equal.  This  ensures  continuous 
dependence  of  the  solution  on  u^  and  uR. 


Figure  1:  Integral  Curves 
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Figure  4A:  Solution  Diagram  for  uL  «  0 

Figure  4B:  Solution  Diagram  for  0^  “  0 

Figure  4C:  Solution  Diagram  for  u^  e  Ai 

Figure  4D:  Solution  Diagram  for  =  0t2 

Figure  4E:  Solution  Diagram  for  uL  e  A2 

Figure  4F:  Solution  Diagram  for  0^  =  ©„ 

Figure  4G:  Solution  Diagram  for  uL  €  A j 

Figure  4H:  Solution  Diagram  for  9^  =  9^ 

Figure  41:  Solution  Diagram  for  u^  e  A4 

Figure  4J:  Solution  Diagram  for  *  -» 
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